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We determine thermal entanglement in mean field clusters of N spin one-half particles 
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.^ ' interacting via the anisotropic Heisenberg interaction, with and without external magnetic 

^^ . field. For the xxx cluster in the absence of magnetic field we prove that only the N = 2 

ferromagnetic cluster shows entanglement. An external magnetic field B can only entangle 
xxx anti-ferromagnetic clusters in certain regions of the B — T plane. On the other hand, 
the xxz clusters of size N > 2 are entangled only when the interaction is ferromagnetic. 
Detailed dependence of the entanglement on various parameters is investigated in each 
^ ' case. 
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1 Introduction 

It is usually stated that at any finite non-zero temperature, thermal fluctuations suppress 
quantum fluctuations and for that reason the latter can exist only at absolute zero or very low 
temperatures. With the progress in quantifying entanglement [T] or quantum correlations, we 
have now a precise method for determining exactly at what temperatures quantum correlations 
cease to exist. This problem has been investigated extensively under the name "Thermal En- 
tanglement" in many body quantum systems 013111112110111112110111111111111211131111113112!. 
The term is usually used to specify the amount of entanglement or genuine quantum correlation 
which exists between two parts of a system when the whole system is in a state of thermal equi- 
librium. One is interested in how this quantity depends on the coupling constants of the system 
and the temperature. Usually there is a threshold temperature above which entanglement van- 
ishes. In some rare cases it is even observed that raising the temperature first increases the 
quantum correlation and then begins to diminish it jH] . In this way by a quantitative treatment 
of the problem we are able to sharpen our intuitive notions on the effect of temperature on 
quantum correlations. 

In addition to these purely theoretical considerations, these works are also motivated by the 
proposals in which spins in solid state systems play the role of qubits jTTl UHl UHl 1201 1211 1221 I23j 
and the fact that thermal entanglement seems to be stable against de-coherence and needs no 
controlled switching for its generation. 

In this regard in the past few years many types of spin systems have been studied, which 
include, the Heisenberg rings with ferromagnetic and anti-ferromagnetic coupling in the ab- 
sence and presence of magnetic fields, small clusters of spins with different types of isotropic 
and anisotropic coupling in homogeneous and in-homogeneous magnetic fields. These are only 
a small sample of the works which have been done in this direction. There are also other works 
which take into account other forms of interactions between the spins. 

Obtaining exact results on thermal entanglement in spin systems on arbitrary lattices is 
extremely difficult, since this generally requires a knowledge of different correlation functions 
which in turn requires a knowledge of the full spectrum of the system. However in certain 
systems with extra symmetries one can obtain some partial and interesting results [Hj. 

In this paper we want to study in some considerable detail thermal entanglement of spins 
in mean field clusters of arbitrary sizes. These are clusters in which every spin is coupled to 
every other spin with equal strength. Such an study has a two- fold motivation. First, mean 
field systems are amenable to exact analytical treatment since as we will see, one can determine 
the entanglement by the sole knowledge of energy eigenvalues, i.e. the partition function, and 
second they are the first approximation that one uses to study other models. Therefore the 
results which are obtained for entanglement in these systems, may be good approximations 
for their corresponding results on other lattices specially when the coordination number of the 



latter is high. 

We would like to stress that in all these considerations we are concerned with the thermal 
state of the system which at zero temperature approaches an equal mixture of various ground 
states, if there is a degeneracy in the spectrum. For the entanglement properties of the ground 
states of xy mean field clusters the reader is referred to PH 123 123 i2Hl 1211 IHO] • 
In these works a general class of mean field models under the name collective models have been 
considered and the entanglement properties of their ground states with particular emphasis to 
its relation with quantum phase transitions have been studied. 

The structure and the results of this paper are as follows: In section |21 we introduce the 
model which consists of a system of A^ spin one-half particles interacting with each other 
and placed in an external magnetic field. The interaction may be of ferromagnetic or anti- 
ferromagnetic type and there may also be present some degree of anisotropy in the couplings of 
spins in different directions. We obtain general formulas for the thermal entanglement of this 
model derived from its partition function. In the following sections we study in detail special 
cases of the model, that is, in section El we consider the isotropic (or xxx model) without 
magnetic field for both ferromagnetic and anti-ferromagnetic interactions, where we show that 
only the anti-ferromagnetic cluster of A^ = 2 spins shows entanglement. In section JH we study 
the effect of magnetic field in the xxx model and show that the ferromagnetic clusters show 
no entanglement while anti-ferromagnetic clusters have entanglement in certain regions of the 
B — T plane, the plane of magnetic field and temperature. Finally in section El we consider the 
anisotropic {xxz model) without magnetic field. For an A^ = 2 cluster both the ferromagnetic 
and anti-ferromagnetic interactions produce entanglement, however remarkably for A^ > 2 we 
see that only ferromagnetic clusters show entanglement and no entanglement develops in anti- 
ferromagnetic clusters. We determine the regions in the A — T plane, the plane of the anisotropy 
and temperature where entanglement is non zero. In all cases we determine the dependence 
of the entanglement in terms of the size of the cluster. Needless to say in a mean field cluster 
the entanglement between any two spins is very low compared to the chain, due to the large 
number of neighbors of a site. Therefore we use the re-scaled concurrence [l^, the ordinary 
concurrence scaled by the number of neighbors to measure the degree of entanglement of a 
cluster. 

2 Interacting spins on mean field graphs 

Consider a cluster of spin ^ particles in an external magnetic field and interacting with each 
other. We take the Hamiltonian as 



J 



N N 

H = J^^ Y. i^i ■ ^3 + ^SizSjz) + BY,Siz, (1) 



where Sia = \<Jia and cTja's are Pauli operators. The couphng constant is taken as jf^ to 
ensure an extensive total energy. We also set J = —1 for ferromagnetic and J = +1 for anti- 
ferromagnetic clusters and without loss of generality B is taken positive. The parameter A 
determines the anisotropy of the interaction. For A = we have an isotropic {xxx system 
which in the absence of external magnetic field displays full SU{2) symmetry. This system is 
exactly solvable since one can rewrite it in terms of the total spin operators Sa '■= J2i Sia in the 
form 



H 



J 



N- 1 



{S ■ S + ASl) + BS, 



(2) 



where S"^ = S"^ + Sy + S^ and we have ignored an additive constant. Note that in studying the 
ground state properties where one can restrict to a fixed spin sector (i.e. S'l + S'^ + S'l = S"^), the 



model is equivalent to the Lipkin-Meshkov-Glick model [26] which has been extensively studied 
in j27| I2HI 121] • However for the thermal entanglement properties of the model one should take 
into account all spin sectors and the above equivalence no longer holds. 

For any system of spins having the symmetry [S^, H] = 0, it has been shown ^J E] that the 
reduced density matrix between two spins is given by 
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where 



u 
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±// + G^ 



^xx ~r ^yy ~r iiy-Jxy ^yx)^ 



(4) 
(5) 



and we have ignored the labels indicating the two sites in question which for the mean field 
cluster can be any two sites. Here Gap and \x are the spin-spin correlation functions and the 
magnetization per site respectively. 



Gap — {SiaSjli), 



(6) 



(Note that there is slight difference between our notations and that of ^ who take Gab 
(o'aO'b))- For such a density matrix the concurrence is given by [Hin] 



C = 2 max 



0,1 z 



-VM+tJ 



from which one can obtain the entanglement of formation (Eof ) from 

Eofip) = -Alog A - (1 - A) log(l - A), 



(7) 



(8) 



with A = i-v^-*^ . We now note that z is real. This can be proved by resorting to the reahty 
of the Hamiltonian and hence the reahty of the total and the reduced density matrices. It is 
worth mentioning that when the interaction strengths between the x, and the y terms are not 
equal, the symmetry \H^ S^ = will no longer hold and the reduced density matrix will not be 
of the form Q. For these cases one can resort to a recent conjecture of ^7\ who provide the 
closed form of the concurrence for general models having permutation symmetry. 



Thus the concurrence can be expressed solely in terms of the correlation functions 
C = 2 max 0, 1 G^x + G 



yy 



G,-.Y - /i2 



(9) 



The magnetization per site /x, and the energy per site e := -^ are obtained from the partition 
function as follows: 



/^ 



1 dlnZ 

'NP dB ' 
1 d\nZ 



N dp 
The correlation function Gzz is obtained from either of the following relations: 
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The other correlation function Gxx + Gyy is obtained from 



yy 
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j{GxX + Gyy + (1 + A)^^ 



liB, 



(12) 



(13) 



where we have used the exchange symmetry between all the pairs in the cluster. Thus from (??) 
and (fTSjl we can determine all the correlation functions which are necessary for the calculation 
of the concurrence. 



The partition function is found from the following formula 



2 O 



Z= Y. Y. e"'^i^(^(^+^)+^'")+^"^<7(iV,5) (14) 

5=0, 1 ra=-S 

where we have ignored a multiplicative constant and the lower limit of the sum begins from 
S = 0{S = |) for even (odd) A^ and 

2S + 1 



^^^'^)=(f-^;fT^TT 



(15) 



is the number of times a spin S representation appears in the decomposition of the tensor 
product of A^ copies of spin | representations. Note that we have to make a distinction between 
odd and even number of sites A^, since for even A^ (odd A^) only integer (half-integer) spins 
appear in the decomposition. We have now set up the required equations and are ready to 
analyze various situations. First we consider the isotropic model in zero magnetic field. 

3 The XXX model in zero magnetic field 

In this case we have B = A = 0. For a finite cluster we have no spontaneous magnetization 
and hence fi = 0. Moreover by rotational symmetry we have G^x = Gyy = G^z = fj- Therefore 
we find from Q that 



C = 2 max 



0, 
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(16) 



where we have abbreviated j to x. We now note that any state with total spin S has an energy 
given by 

Es = j^{S{S + 1)-\n), (17) 

from which we can derive a bound for the variable x = j. This bound is obtained by considering 
the maximum and the minimum values of 5", respectively given by ^ ^ind 0(|) for even (odd) 

A^. 

3 1 

<x < - 



^ A^ = 


even 




^ A^ = 


odd. 


(18) 


ion /(x) := 


- ^ X - 
3 -^ 


- 1 + f satisfies the 



4(A^-1) - - 4 

3 1 
< X < - 

4A^ - - 4 

On the other hand it is readily seen that the function /(x) : 
following bound 

< /(x) only if x < — . (19) 

Combination of (jl8|) and (|19p shows that there is no entanglement for any isotropic cluster for 
A^ > 3 . We can also exclude the ferromagnetic N = 2 cluster by noting that for such a cluster 
e < (sine at infinitely high temperatures e oc tr{H) = 0, consequently at lower temperatures 
e < 0) and J < making the fraction x a positive number. Thus we arrive at the result that 

The only isotropic mean field cluster which can be entangled at non-zero temperature is the 
N = 2 anti-ferromagnetic cluster. 

It had already been shown [H] that in the class of anti-ferromagnetic rings, the A^ = 3 case 
does not show entanglement, here we see that this is a special case of a more general result 
which holds for all mean field clusters, including A^ = 3 clusters as a special case. 




Figure 1: Color online. The re-scaled concurrence for an A^ = 3 xxx cluster, a triangle of 
anti-ferromagnetic spins, as a function of temperature and magnetic field. In this and all the 
other figures, the units are so chosen that B and T become dimensionless. The concurrence is 
a dimensionless quantity. 



4 The xxx model in magnetic field 

In the presence of magnetic field, the model will no longer have the su{2) symmetry and the 
relation G^x = Gyy = Gzz = ^ ^^o longer holds. However one can use (fT^ and calculate the 
partition function which now takes the form 



N_ 
2 



z 



E 

5=0i 



^ ^^,sis+^)'^I^^Mg±jlg^N,S). 



sinhf 



(20) 



All the required correlation functions can be obtained exactly, but can not be written in terms 
of short expressions. The concurrence is obtained numerically. The result is that no thermal 
entanglement develops in the ferromagnetic clusters but for anti-ferromagnetic clusters, there 
is entanglement which generally but not always diminishes by increasing temperature. Figures 
(^ and show the thermal entanglement of A^ = 3 and A^ = 4 clusters as functions of 
temperature and magnetic field. In both of them we see that there are regions in the B — T 
plane where an increase of temperature first increases the entanglement and then tends to 
decease the entanglement. The maximum entanglement exists at zero temperature and inside 
a certain interval of magnetic field values. 

In order to depict the typical information contained in these two figures for different cluster 
sizes, we refer the reader to figures (jHI) and (jlj. Figure (jH} shows the entanglement at very low 




Figure 2: Color online. The re-scaled concurrence for an A^ = 4 xxx cluster, a tetrahedron of 
anti-ferromagnetic spins, as a function of temperature and magnetic field. 



temperatures (T=0.01) for xxx clusters as a function of magnetic field S, for cluster sizes. It is 
seen that by increasing the size of clusters, the interval for which there is entanglement shrinks 
to a vanishingly small interval centered around B = 1. Figure (0} shows the region in the B — T 
plane in which the cluster is entangled. The region becomes smaller and smaller as the size of 
the cluster increases. The right-tilted shape of these regions indicates that in certain intervals 
of magnetic fields one can generate entanglement simply by increasing the temperature. The 
physical reason is that when the magnetic field is high, the ground state is a state with all 
the spins aligned in the direction of magnetic field and hence there is no entanglement at zero 
temperature. Slightly increasing the temperature mixes the entangled excited states with the 
ground state and generates entanglement. 

Finally figure (0) shows the maximum re-scaled concurrence for clusters of different sizes. 
By re-scaled concurrence we mean the concurrence between two sites multiplied by A^ — 1 which 
is the number of neighbors of a given site. It is seen that there is almost no entanglement for 
clusters of size larger than 23. 

This concludes our investigation of the xxx clusters in magnetic field. We now turn to 
another type of anisotropy, namely the xxz clusters in which the anisotropy is not brought 
about by an external magnetic field. 




Figure 3: Color online. The re-scaled concurrence for different mean clusters at very low 
temperatures (T=0.01) as a function of magnetic field. The curves correspond (from top to 
bottom) to (A^ = 2) to (A^ = 12). The units are so chosen that B and T become dimensionless. 
The concurrence is a dimensionless quantity. 

5 The xxz model in zero magnetic field 

For these clusters we find a surprising results: Except for N = 2, the anti- ferromagnetic clusters 
show no entanglement under any condition, regardless of the value A and T. However ferro- 
magnetic clusters show entanglement for all negative values of the anisotropy parameter A. 
We do not report the result for the N = 2 ferromagnetic cluster, since this has been reported 
elsewhere jTH] and focus instead on clusters of arbitrary sizes. 

The entanglement of a ferromagnetic cluster of size A^ = 20 is shown in figure ©. This figure 
is typical, that is all the ferromagnetic clusters have a similar behavior. 



In order to depict the relevant information for clusters of different sizes we refer the reader to 
figures through (fTT|l . Figure © and (fTTHl show the entanglement at very low temperature, 
i.e. T = 0.01, for even and odd size clusters as a function of A. Note the difference in trends, 
that is, for even size clusters as A^ increases the re-scaled concurrence decreases while for odd 
size clusters, as A^ increases the re-scaled concurrence increases. The two figures suggest that 
in the limit A^ — > oo, the re-scaled concurrence approaches a curve which we redraw in fig- 
ure (lllj). The insets in these figures i.e. figures and (|10|) show the regions in the A — T 
plane below which thermal entanglement exists for different clusters. For both even and odd 
size clusters, as the size of the cluster increases, the regions decrease in size but do not shrink 
completely and approach a limiting region. 
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Figure 4: The regions of entanglement in the B-T plane for mean field clusters of size N=3 to 
size N=8 . 



We have to understand two important characteristics of the entanglement of these clusters. 
First why in the limit A — > oo, where the interaction approaches an Ising-like interaction we 
still have entanglement at very low temperatures and second why in the anti-ferromagnetic case 
there is no entanglement at all for any value of A. To this order we present a simple examples, 
namely the = 3 cluster. 

5.1 The N=3 cluster 

For such a cluster, the hamiltonian is 

H = J{si ■ S2 + S2 ■ S3 + Si ■ S3) + JA{sizS2z + Si^Ss^ + S2zS3z) = JS ■ S + JAS^, (21) 

where in the second form, we have ignored an additive constant. Any state of the form \s,m) 
where s and m are respectively the total and the z component of the total spin is an eigenstate 
of energy with energy given by Eg^m = J{s{s + 1) + Am^). The eigenstates are as follows: 
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Figure 5: Maximum re-scaled concurrence for isotropic mean field clusters in magnetic field as 
a function of their size A^. For A^ > 23, the concurrence vanishes. 
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The energies of the states are one of the values ^(15 + 9A), -1(15 + A) and ^(3 + A), de- 
pending on the values of the quantum numbers s and m. 

Figures (jH} and ((Tj) show the spectrum in the ferromagnetic and anti-ferromagnetic cases sep- 
arately. In each case the nature of spectrum depends on the value of A. 
Let us consider the ferromagnetic and anti-ferromagnetic cases separately. 
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Figure 6: Color online. Re-scaled concurrence for an A^ 
function of temperature and anisotropy. 



20 xxz mean field cluster as a 



Anti-ferromagnetic interaction: J = 1 

Here we have the ground state energy Egs as 

i(15 + 9A) for A<^, 

i(3 + A) for ^<A 



E. 



9S 



(31) 



Consequently ioi A < ^, the thermal state at zero temperature is equal to 

p(r = o) = ^(|+,+,+)(+,+,+| + |-,- 



-)(- 



which is clearly a separable state. For ^ < A, the thermal state is an equal mixture of the 
four states with s = ^. A simple calculation gives in this case, the explicit form of the two 
particle density matrix as 
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(32) 



with eigenvalues |, |, | and |. Since in this case we have ^/pp := J p{o-y ® (jy)p*{cry ay) 
p, we find the concurrence of this state to be 
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Figure 7: The spectrum of the A^ = 3 xxz anti-ferromagnetic cluster. Note that I2, 2) stands 
for two doublets. 
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Figure 8: The spectrum of the A^ = 3 xxz ferromagnetic cluster. Note that 13)2) stands for 
two doublets. 



C = max (0, 3-) = 0. 



(33) 



Therefore we have shown that this anti-ferromagnetic cluster is not entangled for any value of 
the anisotropy. 

Here we have used the original form of the concurrence [Jj valid for any density matrix of two 
qubits which states that 

C = max (0, Ai — A2 — A3 — A4), (34) 

where Aj's are the positive square roots of the eigenvalues of pp. 
The ferromagnetic interaction: J = — 1 
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Figure 9: Re-scaled concurrence for anisotropic mean field clusters at nearly zero temperature 
(T=0.01), for clusters of even size, from N=2 to N=20. The inset shows the regions in the A — T 
plane below which there is entanglement. The units are so chosen that T is dimensionless. The 
parameter A and the concurrence C^ are dimensionless quantities. 



In this case the ground state energy Egg is 



E, 



gs 



"^(15 + 9A) 



(35) 



for A < 0, 

, ,.. , .^, for < A 

Consequently for A > 0, the thermal state at zero temperature p{T = 0) is equal to 

p(T = o) = ^(|+, +,+)(+, +,+! + !-, -,-)(-,-, -I) 

which is clearly a separable state. On the other hand For A < 0, the thermal state is an equal 
mixture of the two states with ||, |) and ||, ^). A simple calculation gives in this case, the 
explicit form of the two particle density matrix as 
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Figure 10: Re-scaled concurrence for anisotropic mean field clusters at nearly zero temperature 
(T=0.01), for clusters of odd size, from N=3 to N=21. The inset shows the regions in the A — T 
plane below which there is entanglement.) 

with eigenvalues |, |, | and 0. This gives the concurrence to be 

4 11 
C = max (0,--2-) = -, (37) 

Therefore this cluster is entangled only for negative values of A and for ferromagnetic interac- 
tion. 

Let us see why at zero or very low temperature a ferromagnetic cluster shows entanglement 
in the limit of very large and negative A, despite the Ising like appearance of the interaction. 
The reason is that in this case, the Hamiltonian approaches H = —{S"^ + AS"^) ~ —AS"^. 
Since A < 0, the lowest energy states are those with Sz = ±|, and these states, as shown in 
figure (jHl) are the states ||, ±i) which are highly entangled. One may now ask why the same 
phenomena does not happen in the anti-ferromagnetic case for very large positive A, for which 
the Hamiltonian takes a similar form? The reason is that as shown in figure (|7j) the ground 
state is now a mixture of four S^ = ±| states, and this mixture does not have any entanglement, 
although the states themselves may be entangled. 
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Figure 11: The re-scaled concurrence for an xxz mean field cluster in the limit of large A^ as a 
function of A, at very low temperature T = 0.01 

6 Discussion 

We have done a rather detailed study of the pairwise entanglement properties of a mean field 
cluster of spin one-half particles interacting via the Heisenberg xxz interaction. We have 
considered the ferromagnetic and anti-ferromagnetic interactions and the isotropic {xxx) and 
anisotropic [xxz) cases with and without magnetic field. In each case we have determined 
the dependence of entanglement on various parameters like the external magnetic field, the 
anisotropy parameter, the cluster size as well as temperature. 

It has already been shown that, although the two particle concurrence of the ground state 
vanishes in the thermodynamic limit, the re-scaled concurrence does not vanish and shows 
quite interesting and nontrivial behavior j27l |2Hl 123 EOl • Some of our results (figures Q and 
(jlOp ) shows that as the size of clusters approach infinity, the re-scaled concurrence at non-zero 
temperature approaches a specific function of the control parameters. In a future work we 
will focus on the thermodynamic limit of these clusters to investigate these functions where 
we should also take into account the possibility of spontaneous magnetization and symmetry 
breaking. 
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